Charge transfer through molecular junctions within Redfield theory: 

subtleties and pitfalls 
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Charge transfer through nanoscale junctions connecting metallic leads with quantum dots or 
single molecules is often described within an open system formulation in terms of Redfield theory. 
Under non-equilibrium conditions, the usually invoked rotating wave approximation is not justified 
which may lead to unphysical steady state solutions with e.g. negative populations. In this work 
we explore subtleties and constraints of the approach and thus clarify its applicability in numerical 
calculations. General findings are illustrated for an analytically solvable case of a molecule with two 
electronic states. 
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I. INTRODUCTION 

Electronic transport through nanoscale junctions has 
been in the focus of intensive research for several years 
now, particularly in the context of tailored quantum dot 
structures and in the context of molecular electronics 1 - - — . 
Theoretical descriptions are challenging as they have to 
capture the complex non-equilibrium many-body physics 
due to the interaction of electronic and phonon states 
on the dot/molecule with those in the leads. Several 
strategies have been developed in the past, each with 
its strengths and each with its limitations. Ab ini- 
tio approaches start from the chemical microstructure 
of the dot/molecule- leads compound but have difficul- 
ties to include e.g. strong correlation effects, strong non- 
equilibrium, and inelastic tunneling processes. Efficient 
formulations to describe these many-body phenomena 
work with model Hamiltonians which have been very 
successfully applied for decades in mesoscopic physics. 
Here, theoretical studies of charge transfer have been 
put forward either within the well-established framework 
of single-particle non-equilibrium Greens functions 4 ^— or 
within the context of open quantum systems?- 1 - . The 
former ones consider the dot/molecule as a tunneling 
barrier with energy-dependent transmission, while for 
the latter ones the leads constitute fermionic heat baths 
which must be properly eliminated to arrive at a reduced 
description of the dot /molecule degrees of freedom 1 -^—. 

In both cases, approximate treatments are required to 
access the charge transfer process, thus concentrating on 
various domains in parameter space such as, for exam- 
ple, the single charge sector for weak tunnel couplings 
including inelastic processes, on many body states on 
the intermediate aggregate, or on elastic scattering in 
the coherent regime. However, whenever the fundamen- 
tal physics changes by tuning external parameters (bias 
voltage, gate voltage, coupling, temperature etc.), one 
typically also has to switch to a complementary approach 
to avoid inconsistencies or even unphysical results 1 ^. 

In the context of a reduced description, a powerful tool 
are master equations (ME) based on a perturbative treat- 



ment of the interaction between the relevant agent and 
its surrounding. In case of bosonic heat baths this leads 
to the famous Redfield equation 8 *! 8 , or, with additional 
constraints imposed, to Lindblad theor y 11 : 19 . Things be- 
come more intricate though when the environment drives 
the agent out of equilibrium and towards a steady state as 
in the case of voltage-biased tunnel junctions. Namely, 
under these circumstances the usually invoked rotating 
wave approximation (RWA) which decouples the dynam- 
ics of diagonal elements of the reduced density (popula- 
tions) from those of the off-diagonal ones (coherences) is 
not justified. Mathematically, the existence of complete 
positive densities is then no longer guaranteed. Apart 
from the failure of this formal concept, the question arises 
to what extent or under which conditions physical accept- 
able steady states, i.e. those with positive populations, 
are actually supported by MEs in practical applications. 

To shed some light on this issue, is the intention of 
the present work. Some of these and related questions 
have been addressed already previousl y 12 : 14 ' 16 , but to 
our knowledge no consistent deeper analysis has been 
given yet. As we are interested in generic properties, 
we consider a strongly simplified set-up compared to re- 
alistic junctions, namely, a single dot/molecule with dis- 
crete electronic states interacting with broadband metal- 
lic leads (Sec. |TT]). Integrating out the lead degrees of free- 
dom brings us to the conventional Redfield equation with 
the Redfield tensor as its central ingredient (Sec. IIIip . In 
steady state the underlying Markov approximation al- 
ways applies so that the proper expansion parameter for 
the Born series is not immediately obvious. We formu- 
late constraints that the Redfield dissipator has to obey 
for physical steady states to exist as stationary solutions 
of the Redfield equations. It turns out that already the 
zero- voltage situation must be treated with care since the 
density operator of the compound may significantly differ 
from the product state of its bare constituents. General 
findings are illustrated by considering a molecule with 
only two electronic levels which allows for a complete an- 
alytical solution, Sec. IIVI In the linear response regime 
contact is made with the non-equilibrium Greens func- 
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tion formulation (Sec. |V|) . Conclusions summarize our 
results and specify under which conditions a Redfield ap- 
proach is applicable for charge transfer (Sec. I VI)) . 



II. SYSTEM AND REDFIELD FORMULATION 

We consider a minimal model consisting of a quantum 
system with Hamiltonian 



]e a \a)(a\ 



(1) 



which is positioned between two fermionic leads I = L, R 
with 



(2) 



where the sum is over a multi-index a = (k,l). In the 
sequel, we will keep this notation and use Greek letters 
for lead states and Latin letters for dot states. As usual, 
electronic states in the leads are Fermi distributed 

(40 = z — < \ i i = Z( e ") 



exp(/3e a ) 
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at inverse temperature /? and chemical potential Vi so 
that e a — efc — Vi with the single particle energy = 
h 2 k 2 /2m. 

The tunneling of electrons between leads and dot is 
described by 



(3) 



Here, the operator $ a removes an electron from the dot 
to state a. 

At vanishing bias voltage the full equilibrium den- 
sity operator can be written as W — exp[— (3(H + 
Hb + Hi)]/Z, where the partition function is given by 
Z = Tr{exp[- /3(H + H B + Hi)]}. The reduced den- 
sity operator of the dot is given by tracing out the 
leads, i.e., p = Tib W. In the regime of weak lead-dot 
coupling, one could then assume that to leading order 
W fa Wp ro d = po ® W B ~ e _/3ff ° e~^ Hs . Such an ansatz 
is, of course, only justified as long as higher order terms 
are sufficiently small. However, this is not always guar- 
anteed. Namely, exploiting the Baker-Hausdorff formula 
one has for the next order contributions (see Ref. [22] ) 



W prod • exp(-f3Hi) ■ exp ( ^-[H Q + H B , Hi 



's V (^[[H + H b ,Hi},Hi 



■exvl^[[Hi,H + H B ],H +H B ]\. (4) 

Due to the fact that Hb oc c?c and Hj oc T (c + c>) 
with typical coupling parameter T, the double commuta- 
tor [[H B , Hi], Hi] is proportional to T 2 times the iden- 
tity in the lead-subspace. Taking the trace may thus 



yield an ill-defined reduced dot operator depending on 
the operators $ a but no matter how small the coupling 
is. The physical reason for that is easy to understand: 
Hybridization between lead and dot states immediately 
occurs for T ^ such that for e.g. Fermi levels in the 
leads lying above the lowest unoccupied level of the dot 
(LUMO), the ground state of the reduced lead-dot com- 
pound is orthogonal to that of the bare dot. Accordingly, 
microscopic approaches such as density functional theory 
(DFT) calculate the equilibrium state of the coupled dot 
or molecule by including part of the leads (the so-called 
super-molecule )2. 

At a finite bias voltage, charge transfer sets in and we 
are interested in the corresponding steady state. This is 
approached during the time evolution of the compound 
starting from an initial state, e.g. an uncoupled equi- 
librium state of system and leads. One way to do so, 
is to derive an equation of motion for the reduced den- 
sity of the dot for weak tunnel couplings. Here "weak" 
is often simply understood as a small coupling param- 
eter r. One purpose of this work is to reveal more 
precisely in which sense such a perturbative treatment 
actually applies. We thus proceed with a conventional 
Born-Markov approximatio n 7 ' 8 ' 10 and start by expand- 
ing the Liouville-von Neumann equation of the full sys- 
tem —iK^xf — [W, H] up to second order in Hi. As a 
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result one obtains the well-known equations' 

f = ~[p,H ]- ± j™ dTTr B {[Hi,[Hi(-r),W(t)]]} 

= ~[p(t),H ]+Rp(t), (5) 

where the Redfield tensor R accounts for the impact 
of the interaction Hamiltonian Hi. The time depen- 
dence in the integral follows from the interaction picture 
Hi(-t) = UqHjUI = exp[-f{H +H B )]Hiexp[f(H + 
H B )]- The central assumption in ([5J are macroscopic 
reservoirs such that W(t) = p(t) (g> exp(— /3H B ). 

Steady states are now inferred as the stationary so- 
lutions to the ME. The Markov approximation in ([5]), 
i.e. a sufficiently short-ranged reservoir correlation func- 
tion, is exact in this case. A subtlety arises, however, 
namely, whether one searches for solutions in the inter- 
action or in the Schrodinger picture^. The first one is 
determined by Rp s t = 0, while the second one follows 
from Rps.st — ~ j^[ps,st, Hq]. It is well-known though 
that the latter ones never carry any coherences, i.e., fi- 
nite off-diagonal elements of ps.st, in a consistent second 
order perturbation theory and as long as the energy spac- 
ings of the states \a) are larger than the tunnel coupling 
constant T = J2k^klt ( see H)- ^ ne argument is 
simple: During the time evolution off-diagonal elements 
(a\ps(t)\b) oscillate with frequency uj a b — (e a — €{,)/% and 
for large times are thus washed out. This is not the case 
in the interaction picture where effectively one works in 
a rotating frame and may thus be able to keep reservoir 
induced coherences between dot/molecule states. 
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Once the steady state is at hand, the most interesting 
quantity one can calculate is the current through lead I: 

(It) = e lim (N t ) = e lim Tr{NiW(t)} 

t—>oo t—t-oo 

= -jTtiNilHo,/*]} 
e f°° 

-T2 dTTriNtlHi^Hrt-T^W*]]} 
ft Jo 

= j drTr{p st ([J ()J ff J (-r)]) B } (6) 

with the lead I number operator Ni = J2 a ^a,l c a c a (^a,l is 
the Kronecker symbol which is only non-vanishing if the 
lead index in a is equal to /) and the current operator 

Jt = [Ni , Hi] = E a <W7a( $ a4 - $ l c ")- In turns out 

that in state representation the current is determined by 
coherences of the steady state (see below). 

Unfortunately, the steady state operator p st obtained 
from ([5]) does not have to be positive though, as we show 
now. For this purpose, we consider the fermion propaga- 
tors of the leads 



ia (r) := (4(r)c Q (0)) B = e^ f{e a - Vi) (7) 



and 



: = ( c L(- T ) c Q (0))s = exp(--e Q r)/(e Q - Vi), 

k q (t) := (c a (T)cl(0)) B = exp(-^e Q r)/(Vi - e Q ), 

n 

<(r) ~ (c a (-r)ci(0)) B =exp(^e a r)f(V l -e a ).(8) 

Now, after cyclic rearrangement of the c a and c* a under 
the trace, one arrives at 

f°° - 

Rp = dry"^[Ka*a*L(-r)p + ««* + $a(-r)p 

Jo a ft 

+ <p*a(-T)#t + R* aP & a {-T)$ a - K*$p$ a (-r) 

- <$ q( 0$ q (-t) - K a $ a (-7>$+ - K a $t (-<r)p$ a ] , 

(9) 

where in state representation 3?a,ai>( — t) = 
(a\$ a (-T)\b) = $ Q!a& exp(ia; Q br) with huj ab = 
(e a - e b )/h. 

Now, according to Lindblad mathematically complete 
positivity of p s t is only guaranteed if the dissipator R is 
of the for m 11 ' 19 



In the above ©, one then must have 

= L^L. (10) 



This can be archived if (i) one applies in addition a sec- 
ular approximation 8 - (also called rotating wave approx- 
imation, RWA) so that the exp(— ie a r) in the k can- 
cels the exp(— ifwj a bT) in the $(— r), or if (ii) one drops 
the r-integration (this leads to the singular coupling 
approximation 8 ) or if (hi) one 'symmetrizes' the range 
of integration dr — > (1/2) J^^dr. However, under 
non-equilibrium conditions (i) may not be justified (see 
e.gi2£), (ii) is in general not realized, and (hi) is an ad 
hoc procedure. This is in contrast for e.g. a bosonic heat 
bath interacting with an autonomous system, where the 
asymptotic state of the reduced system is the bare ther- 
mal equilibrium ~ e _ ^ ff ° so that coherences (off-diagonal 
elements) die out during the time evolution. The RWA or 
related approximations are then applicable under much 
milder conditions. Here, if none of these reductions work, 
one may wonder if (or under which conditions) at least 
positivity (all populations (a|p st |a) > 0) survives so that 
stationary solutions of © yield physical steady states. 
This will be the central subject of the remainder of this 
work. 



III. THE REDFIELD DISSIPATOR 

In the following, in state representation we think of 
R = R(ab,cd) in © as a matrix with 'columns' (ab) and 
'rows' (cd) which then can be written in the form 

c' ,a 

+^c]T(^ c ,aW 6 £- t , d + $a )dc '< c , 6 ££ c , rf ) 
c' .a 

+$ a , db $l ac (e£ m + £™ c )] . (11) 

Here we introduced the self energy 



S a ab = la I dr K a (-r) exp(iuj ab T) (12) 
Jo 

which is identical to the Keldysh inscattering func- 
tion [2l| if the transition from state b to a moves 
an electron from lead to dot, respectively. The 
corresponding outscattering function reads S° ut = 
7 a / °° dTK a (— t) ex.p(iuj a bT). Thus, the time evolution 
of the populations p aa is governed by 



Raa.cd = S ad £ ^^cE*^ + £ ^ a , dc ^ a ,c'a^7d 
a,c' a,c' 

" E ®lda®a C (KM + KTc) + n ^ out ) -(13) 



This matrix is in general complex- valued and couples not 
only populations, but also populations and coherences. 
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Now, due norm conservation the row trace of R is zero 
{Property 1), i.e., 

a, a c' 

a, a 

+{in out} = . (14) 

This in turn means that the rank of R is not full so that 
Rpst = always has a non-vanishing solution. 

In particular the stationary current through the dot 
(J6j> is thus given by 

/>oo 

(I,) = -5 / drT>.{pi rt <[J,fl>(-T)]) I } 
Ji Jo 



J2 ^>lArt,a&($a,6c$a,ca S a,6a 



(15) 



A. Self energy 

To further elucidate the properties of the Redfield dis- 
sipator, we explicitly calculate the self-energies S. Ac- 
cordingly, the wide band limit is taken and coupling con- 
stants are assumed to be independent of the fermion state 
7 a = 7; with bands of constant density of states Di in 
the leads. As a consequence, one has 

a 

POO 

= / dTY,Kale* eaT (cU-T)c k (0))ie-^ 

= TiJ dTe^-^y dee-* eT f(e) (16) 
with Ti = Dijf. The second integral yields 



J de e-i-/(e) 



/3sinh^ ' 



(17) 



where r includes a positive imaginary increment i/uj c . 
Therefore, with the abbreviation uji. a b — co a b — Vi/h we 
have 



J Lab 



00 n lU,„ h T 



dr- 



7TT 

hp 



fi>P Ji/u c sinh- 

/OO J^L 
-M n=0 



h/3 

=: E(uj ab -Vi), 



7a + *(__ t _j + _ + I n f_ 



with Matsubara frequencies v n = 2im/Tij3, Euler's con- 
stant 7s, and in the limit where u c by far exceeds in- 
trinsic energy scales of the dot. Physically, Tiuj c corre- 
sponds to a cut-off energy in the leads and must thus 
assumed to be large in order to be consistent with the 
Markov- approximation and the broadband limit. For 
the real part one regains the Fermi distribution SJ n = 
TrTif(hu}^ab)- The imaginary part S m (u>) is symmetric 
in uj and exhibits a logarithmic dependence on the cut-off 
uj c . Due to S°^ t (cj) = £*"(— w) in- and out-self energies 
contain identical imaginary parts E ln . Note that a sym- 
metrization of the r-integration in Q as discussed above, 
would only yield real-valued self-energies. For a bosonic 
heat bath the corresponding imaginary part is known as 
reservoir induced Lamb-shift which differs from the one 
obtained for the fermionic self energy only by the addi- 
tional term of 1/2 in the argument of the ^-function^^. 



B. Steady state 

We now search for a stationary solution Rp s t = for 
the density matrix in interaction representation, which, 
thanks to Property 1 always exists and has a trace of 
1. To be a physical density, p st also has to be positive 
definite. Equivalently, \p a b\ 2 < Paapbb which means that 
all paa > 0. This in turn forces the row vectors of R 
to be orthogonal on the subspace p aa > 0. Accordingly 
(Property 2): Stationary solutions of Rp s t — are not 
physical steady states if 



3a6 Rab.cc 

> 0. 



(19) 



(18) 



In other words, rows of R must contain at least one neg- 
ative element for a physical steady state to exist. 

In the secular approximation, i.e. by dropping all parts 
of R which mix coherences and populations, the existence 
of a physical solution p s t is always guaranteed. Then, the 
conservation of the total probability renders the rows of 
R to be linear dependent so that each column has a sig- 
nature of the form H (the ith column has its 

plus at the ith position) in the subspace of the p QQ , thus 
violating (p~9|) . For the full Redfield dissipator (fTTj) , how- 
ever, "off-diagonal" rows of R do not need to obey a cer- 
tain signature and the existence of a physical acceptable 
solution depends on details of the problem under inves- 
tigation. To analyze this analytically more carefully, we 
consider a two level system in the next section. 



IV. EXAMPLE: TWO LEVEL SYSTEM 

The system consists of two electronic sites |1) and \2) 
with energies e±/2 — ±5/2 and a hopping element A be- 
tween them. Site |1) couples to the left lead with a cou- 
pling strength and \2) couples to the right lead with 
a coupling strength of 7^. Diagonalizing the uncoupled 
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Hamiltonian yields the singly occupied energy eigenstates 



A. Redfield tensor 



-) = cos^|l)+sin^|2) 
-> = - S in^|l)+cos^|2) 



(20) 



Now, due to charge quantization, in the relevant den- 
sity matrix 



with eigenvalues H \±) = ±A|±) where A = \/5 2 + A 2 
and the mixing angle is tan((/>/2) = A/ (A + 5). Hopping 
to and from the dot (with the empty dot denoted with 
|0) at energy eo) is described by the Hamiltonian Q with 
matrix elements 



P = (poo,P++,P — ,P+-,P-+) , 



(22) 



sectors with a different number of charges on the dot do 
not mix. With the shorthand notation 



$l,o+ = cos- , $ fli0 + = sin- 
$L,o- = -sin- , $r,o-=cos- 



and $;, aa = $ /:+ _ = 0. 



$^ := $ +0 $ t + , $ 2 _ := $-o$J_ = , n := $ +0 $J_ , 

(21) (23) 
and using the fact that E±o = £±o an< ^ Eo± = Eq± the 
Redfield matrix (JTTJ) R = R' + iR" reads 



/ 2<^E' +0 + 

-2$|S' +0 
-2$ 2 E' 



2$ 2 _E'_ 



V 



-n(s'_ 
-n(s' +0 



-0 

S'_o) 



-2$5.5y 0+ 

2<^E' 0+ 




nE^ 



0+ 



-2$2_E^_ 



2$ 2 _E£ ) _ 
HE[,_ 
HE[,_ 



n(E^ + E^ + ) -n(E(, + + E^_) \ 

rrn/ rrn/ 



nE' 



0- 



nE' 



$iEf 



o+ 



/ 










n(E£ + - e 


























-HE^ + 


n(E% - 




HE'oV 


-ns' '_ 




\n(E%- 


£%) 


-nEg+ 








n(E£_ 



-ns£_ 
ns£ + 
o 



^0+ 
\ 



■ $ 2 _E[,_ / 



v" 

^0 



(24) 



where the lead index I applies to all entries of the ma- 
trix. In the diagonal of R", the imaginary part of the self 
energy gives rise to a Lamb shift 

^ LS = ]T(<f 2 ,+E2; 0+ - $2_e;; _) , (25) 
1 

which can simply be absorbed into a re-definition of Hq 
and will thus be neglected in the sequel. 

Since the Hamiltonian is hermitian, we have R^ >aa = 

R' aa + iR'a a = R*_ + aa - Hence, together with Hr = 
—Hl = 2 sin^ one finds 

R'aa = l)'(^-0i + ^+0/) — ^O+i' — 

I 

Ka = ^E(- 1 )'( S -o i - S +o^o +i ,-4'- i )(26) 
I 

which may have any signature, especially at non-zero bias 
Vj, = Vl — Vr ^ or for asymmetric couplings to the 



leads Y l ^ Tr. One always obtains an acceptable steady 
state for Vb = 0. This in turn motivates a perturbative 
treatment around the zero-bias situation (see second last 
section). 



B. The case of zero temperature 

The entries in the Redfield tensor simplify considerably 
in the zero temperature limit T — > 0. The self energy 
reduces to 

e£ 6 = 7rr,e(-w Il06 ) 

<™ b = -r^ + hiK^MI) (27) 

with the step function 0(w). Accordingly, there are six 
generic cases with respect to the voltages since we can 
always assume that Vl > Vr and e + > e_, namely, 
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/l-o 


Jl+q 


fn-o 


fn+o 


Pst 














+ 


1 











? 


1 





1 





+ 


1 


1 










1 


1 


1 





? 


1 


1 


1 


1 


+ 



where Y,(Tiujba) 



-'l, be 



v b =o 



The conductivity is of order 



Here, we used for the Fermi functions the abbreviation 
flab = f{^i,ab) = S(-w;,afe)- in the last column, cases 
marked with a '+' always yield acceptable steady states 
(positive populations), cases with a '-' yield unphysical 
steady states (at least one population is negative) and 
cases with a '?' lead to unphysical steady states if Tl ^ 



V. PERTURB ATIVE TREATMENT FOR 
SMALL BIAS 

For Vb — steady states are always acceptable sta- 
tionary solutions. Hence, one may wonder if at least 
a perturbative treatment of the the full Rcdfield ap- 
proach for low voltages always applies. For this pur- 
pose, we expand the Redfield matrix around V> = 0: 
R w Rq + Vb§$- = Rq + V b 5R. The steady state density 
is then given by p s t = po + V>Sp with i?oPo = 0. Note 
that po m general differs from the density of the bare 
dot and thus accounts also for lead-dot hybridization as 
we discussed around (j4|). The first order correction is 
determined from 



Ro Sp + SR po = . 



(28) 



As long as we are not near resonances such that Vb differs 
sufficiently from transition energies TiLo ao , do the columns 
of SR lie in the space spanned by the columns of Rq. 
The solution for 5p is then simply calculated. Close to a 
resonance though, (|2"8)l has no solution meaning that the 
Redfield scheme suffers from exactly the same difficulties 
as the non-equilibrium Green's function formalismii. 

The steady state current (|15|) away from resonances is 
easily obtained. We have (/;) = V b Gi with G = SpT, + 
PoST,. Since G is an observable and po is real, only the 
real parts of 5S contribute to G (the imaginary parts 
are symmetrized out), which, however, vanish away from 
resonances. In addition, as seen from (I28|) . in this regime 
the main contributions to SR come from the imaginary 
part R". Thus, Sp consists almost entirely of coherences 
and the conductivity reads 



out 
bl 



a be 



abc 

-$J c $ ca S(-nw ba )j , 



i R in the lead-dot coupling and thus captures sequen- 
tial charge transfer only but keeps lead induced coher- 
ences between the dot states. We recall that Sp is deter- 
mined from a physical state po via (|28[) and not from the 
density of the uncoupled dot. 



VI. CONCLUSIONS 

We are now in position to summarize the findings of 
the previous sections and to draw conclusions. These 
refer to general situations and may not necessarily apply 
to specific set-ups. 

(i) The Redfield approach provides always physical 
steady states if imaginary parts of the Redfield tensor, i.e. 
imaginary parts of the self-energies {jTHJ , are neglected. 
However, since in general real and imaginary parts are 
on the same order of magnitude, this procedure (even if 
appealing for practical purposes) is not justified, (ii) In 
steady state the Markov approximation always applies 
independent of the couplings r^, T^. (iii) At zero volt- 
age, the steady state cannot be taken as a product state 
consisting of the bare dot and the bare lead densities, 
but must be calculated from the full Redfield equation 
which then always provides physical states, (iv) For fi- 
nite bias voltages, the only consistent use of the Redfield 
formulation for steady states is to apply an additional 
perturbative expansion around zero bias voltage in the 
interaction representation; this procedure basically co- 
incides with a linear response treatment. The correct 
expansion parameter for the formulation is then not just 
given by the coupling constants Tl,Tr. (v) In the high 
temperature limit hj3wi ra b <C 1 with £J; ja & ^ 0, the treat- 
ment is justified if /3(Tl + T R ) -C 1. Even in this regime 
must a finite bandwidth uj c be kept in (I18[) to avoid diver- 
gencies in the imaginary parts of the self-energies, (vi) In 
the low temperature domain, the low- voltage expansion 
applies if 



\hoj ab \ » k B T, V b , \T L - T R \, (T L + T R ) 



(30) 



(29) 



with the net bias voltage Vb = \Vr — Vl\. 
Here, the inequalities result from the conditions for 
V, dZjdVly^ ~ d^{l/2-i^ a b/vi)/dVi\ Vi=0 to 
be sufficiently small. In particular, the Redfield approach 
is thus constraint to voltages away from resonances and 
to weaker couplings with limited asymmetry. Outside 
this domain, numerical solutions of the full Redfield equa- 
tions © lead to unphysical steady states: The numerics 
includes all higher order voltage contributions contained 
in the second order (~ V 2 ) Born-Markov expansion while 
corresponding voltage contributions contained in fourth 
and higher order (~ T n ,n > 4) are neglected. This is 
inconsistent and can then also not be cured by going be- 
yond the conventional T 2 expansion. 

The physical interpretation is that no matter how small 
the coupling to the leads may be, in the asymptotic long 
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time limit where the steady state exists, correlations be- 
tween agent (dot) and its environment (leads) may drive 
the compound far away from its bare structure. This can 
happen either when external energy supplied by a voltage 
source induces resonances or when large asymmetries in 
the agent-environment coupling induce strong quantum 
fluctuations and coherences between energy eigenstates 
of the bare agent. The Redfield approach is only consis- 
tently be applicable in domains of parameter space where 



these processes are suppressed. 



Acknowledgements 

We thank M. Grifoni, C. Timm, M. Thoss for valuable 
discussions. Financial support was provided by the DFG 
through SFB569 and the German-Israeli Foundation. 



1 C. Joachim, C.K. Gimzweski, A. Aviram, Nature 408, 541 
(2000); A. Nitzan, N.A. Ratner, Science 300, 1384 (2003). 

2 G. Cuniberti, G. Fagas and K. Richter, Introducing Molec- 
ular Electronics, (Springer, New York, 2005). 

3 J. C. Cuevas and E. Scheer, Molecular Electronics: An 
Introduction to Theory and Experiment, (World Scientific, 
Singapore, 2010) and references therein. 

4 Y. Meir and N.S. Wingreen, Phys. Rev. Lett 68, 
2512(2002). 

5 M. Galperin, M.A. Ratner, A. Nitzan, J. Phys. Cond. Mat. 
19,103201 (2007). 

6 A. Gogolin, A. Komik, |cond-mat/0207513| (2002); D. 
Mozyrski, M.B. Hastings, I. Martin, Phys. Rev. B 73, 
035104 (2006). 

7 U. Weiss, Quantum dissipative systems, (World Scientific, 
Singapore, 2008). 

8 H. P. Breuer, F. Petruccione, The Theory of Open Quan- 
tum Systems, (Oxford University Press, Oxford, 2002). 

9 S.A. Gurvitz, Phys. Rev. B 57, 6602 (1998). 

10 A. Mitra, I. Aleiner, A. J. Millis, Phys. Rev. B 69, 245302 
(2004). 

11 R.Alicki, K.Lendi, Quantum dynamical semigroups and ap- 
plications, Springer, Lecture Notes in Physics 286 (1987). 



2 C. Timm, Phys. Rev. B 77, 195416 (2008). 

3 A. Donarini, M. Grifoni, K. Richter, Phys. Rev. Lett. 97, 
166801 (2006). 

4 M. Leijnse and M. R. Wegewijs, Phys. Rev. B 78, 235424 
(2008). 

5 C. Timm, Phys. Rev. B 83, 115416 (2011). 

6 S. Koller, M. Grifoni, M. Leijnse, M. R. Wegewijs, Phys. 
Rev. B 82, 235307 (2010). 

7 F. Evers and K. Burke, Pride, Prejudice and Penury of ab 
initio transport calculations for single molecules in CRC 
Handbook on Molecular Nanoelectronics (ed. S. Lyshevski), 
2007. 

8 A. G. Redfield, IBM J. Res. Dev. 1, 19 (1957). 

9 G. Lindblad, Comm. math. Phys. 48, 119 (1976). 

:0 J. P. Pekola, V. Brosco, M. Mottonen, P. Solinas, and A. 
Shnirman, Phys. Rev. Lett. 105, 030401 (2010). 

1 S. Datta, Electronic transport in mesoscopic systems, Cam- 
bridge University Press, 1995. 

;2 J.A. Oteo, J.Math.Phys 32, 419 (1991). 

!3 V. Gramich, P. Solinas, M. Mottonen, J. P. Pekola, J. 
Ankerhold, Phys. Rev. A 84, 052103 (2011). 



